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i  scaling  fEinctions.  The  qiiestion  of  whether  the  Integer  translates  of  the  scaling 
!  function  constructed  from  a  prescribed  scaling  sequency  In  the  st.andard  way  are 
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■:  arguments  which  are  currently  in  the  literature  art?  very  complicated. 
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Report  Summary 


The  principal  investigator,  W.  R.  Madych,  participated  in  the  AMS-SIAM 
Summer  Seminar  on  the  mathematics  of  tomography,  impedance  imaging, 
and  integral  geometry  held  at  Mt.  Holyoke,  June,  1993.  He  gave  a  lecture 
on  his  recent  work  in  wavelets  and  sampling  at  this  seminar. 

In  the  preceding  reports  we  outlined  several  lines  of  investigation  which 
we  are  pursuing.  The  following  work  has  been  completed  since  the  last  report: 


•  The  recovery  of  irregularly  sampled  band-limited  functions  via  tempered 
splines.  In  this  article  we  show  that  band  limited  functions  can  be 
recovered  from  their  values  on  certain  irregularly  distributed  discrete 
sampling  sets  as  the  limits  of  the  piecewise  polynomial  spline  inter- 
polants  when  the  order  of  the  splines  goes  to  infinity.  This  is  significant 
extension  of  the  classical  case  when  the  sampling  set  is  a  lattice  which 
was  considered  by  L.  CoUatz,  W.  Quade,  I.  J.  Schoenberg,  and  others. 

•  Orthogonality  criteria  for  compactly  supported  scaling  functions.  The 
question  of  whether  the  integer  translates  of  the  scaling  function  con¬ 
structed  from  a  prescribed  scaling  sequence  in  the  stand^d  way  £ire 
mutually  orthogonal  is  quite  subtle.  The  various  conditions  and  the 
supporting  arguments  which  are  currently  in  the  literature  are  very 
complicated.  In  this  article  we  give  new  simple  proofs  of  several  crite¬ 
ria  for  the  orthogonality  of  the  integer  translates  of  a  scaling  function. 


Both  articles  are  included  in  this  report. 

At  present  we  are  preparing  several  papers  detailing  our  grant  related 
work.  The  topics  include  (i)  certain  questions  concerning  irregular  sam¬ 
pling  of  signals,  (ii)  results  concerning  the  breakdown  the  so-called  scaling 
functions  into  more  elementary  building  blocks,  emd  (iii)  results  concerning 
multivariate  wavelets  and  tilings  of  i?". 

Finally  we  mention  that  in  addition  to  the  ongoing  work  we  we  have 
completed  twenty  one  articles  on  various  aspects  of  wavelets  and  their  appli-  — 
cations.  All  of  these  have  been  included  together  with  detailed  summaries  in 
the  five  semi-annual  and  annual  technical  reports  which  have  been  submit-  ^ 
ted.  As  of  this  date  fifteen  have  been  published  or  accepted  for  publication  , 
in  various  scholarly  journals  and  books.  :  on. 
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THE  RECOVERY  OF  IRREGULARLY 
SAMPLED  BAND-LIMITED  FUNCTIONS 
VIA  TEMPERED  SPLINES 

Yu.  Lyubarskii  and  W.  R.  Madych 


Abstract 

We  show  that  band  limited  functions  can  be  recovered  from  their 
values  on  certain  irregularly  distributed  discrete  sampling  sets  as  the 
limits  of  the  piecewise  polynomial  spline  interpolants  when  the  order 
of  the  splines  goes  to  infinity.  This  is  an  extension  of  the  classical  case 
when  the  sampling  set  is  a  lattice  which  was  considered  by  L.  CoUatz, 
W.  Quade,  I.  J.  Schoenberg,  and  others. 


1  Introduction 

1.1  Overview 

This  paper  concerns  univariate  splines  of  even  order  and  band-limited  func¬ 
tions.  The  main  result  asserts  that  functions  in  the  classical  Paley- Wiener 
claiss  PW,,  that  is,  those  square  integrable  functions  whose  Fourier  trans¬ 
forms  are  supported  in  the  interval  [— x,  xj,  can  be  recovered  from  their  values 
on  certain  irregularly  spaced  sampling  sequences  via  the  formula 

lim  Sk{x)  =  /(x) 

fc— *00 

where  s*  is  the  spline  of  order  2k  which  interpolates  /  on  such  a  sequence. 
For  more  precise  statements  see  Theorems  4  and  5  in  Section  4.3  below.  The 
sampling  sequences  {x„}  considered  here  are  those  for  which  the  correspond¬ 
ing  exponential  functions  constitute  a  Riesz  basis  for  1^([— x,7r]). 

1.2  Background,  motivation,  and  objectives 

If  /  is  a  function  in  the  Paley- Wiener  clas,s  then  it  is  uniquely  determined  by 
its  values  on  the  integer  lattice  Z  and  can  be  recovered  from  these  values 


I 


•  • 
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via  the  VVhittaker-Kotelnikov-Shannon  sampling  formula 


(1) 


m 


n=— oo 


fin) 


sin  x(t  —  n) 
/r(x  —  n) 


For  example  see  [27]. 

Such  f's  are  aJso  uniquely  determined  and  can  be  recovered  from  their 
values  on  certain  irregular  sampling  sets  {x„}  =  Xq  <  Xi  <  ...}. 

In  particular  if  {a:„}  is  such  that  the  collection  of  functions  iy7n(0  =  n  6 

Z,  is  a  Riesz  basis  for  i^([— t,  jrj)  then  the  natural  substitute  for  (1)  is 


(2) 


/(x)=  ^  /(x„) 


S(x) 


S'(x„)(x  -  x„) 


where  S(z)  is  the  unique  entire  function  of  exponential  type  tt  whose  zero 
set  is  {x„},  see  [19].  For  more  references  see  the  discussion  in  Section  5. 

Because  siniri  as  well  a.s  S(x)  fail  to  belong  to  L^(—oo,oc)  the  cardinal 
functions 

sin7r(x  — n)  S{x) 

"  ^iict  .  .  "  “ 

t(x  -  n)  S'(x„)(x-x„) 

have  relatively  slow  decay  as  x  — >  ±oo.  This  makes  formulas  (1)  and  (2) 
somewhat  unstable.  A  perturbation  at  n  or  x„  will  have  significant  influence 
at  points  i  far  from  n  or  x„  respectively. 

Schoenberg  used  cardinal  splines  to  stabilize  (1)  in  [23].  Subsequently,  he 
showed  that  if  /  is  in  the  Paley-Wiener  class  PW,  and  s*  is  the  piecewise 
polynomial  spline  of  order  2^'  which  interpolates  f  on  Z  then 


(3)  lira  Sfc(x)  =  f(x). 

fc— OO 

uniformly  in  i.  The  spline  enjoys  the  representation 

(4)  Sfc(i)  =  f{n)^k{x:  -  n), 

rt=— x> 


where  Afc(x)  is  the  corresponding  cardinal  spline  function.  This  fundamental 
spline  A*  is  the  unique  spline  of  polynomial  growth  of  order  2k  with  knots 
on  Z  which  satisfies  Afc(n)  =  5o.n,  n  €  where  6o_„  is  the  Kronecker  delta. 
Since  A*  heis  exponential  decay  at  ±oc,  see  [23],  (3)  and  (4)  represent  a  more 
stable  method  of  approximating  /  than  (1).  For  extension  to  wider  classes 
of  band-limited  f's  see  [24.  22,  18).  The  results  make  essential  use  of  the 
lattice  structure  of  the  sampling  set  Z,  particularly  the  Poisson  summation 
formula. 


It  should  be  mentioned  that  the  use  of  splines  in  a  summability  method 
for  the  recovery  of  regularly  sampled  band  limited  functions  goes  at  least  as 
far  back  as  the  work  of  Quade  and  Collatz  [21]  where  a  variant  of  (3)  was 
established  for  certain  trigonometric  polynomials  /.  See  also  [24,  page  103]. 

In  the  case  of  irregularly  spaced  knot  or  sampling  sequences  {x„}  that 
satisfy  certain  natural  conditions  the  corresponding  cardinal  functions  which 
build  the  interpolating  piecewise  polynomial  splines  are  also  known  to 
enjoy  exponential  decay,  see  [3].  Convenient  ways  of  evaluating  such  splines 
Sk{t)  in  terms  of  the  data  {/(x„)}  can  be  found,  for  example,  in  [4]. 

In  view  of  this  situation  it  is  natural  to  ask  whether  the  spline  summability 
method,  which  makes  essential  use  of  the  lattice  structure  of  the  sampling 
set  Z,  can  be  extended  to  the  case  of  irregularly  sampled  data.  As  indicated 
in  the  overview,  Section  1.1,  we  found  the  answer  to  this  query  to  be  yes.  To 
obtain  this  answer  we  applied  the  technique  of  mean  periodic  continuation 
with  respect  to  a  given  basis  of  exponentials  to  the  study  of  splines  with  a 
biinfinite  knot  sequence.  The  goal  of  this  article  is  to  outline  the  theory  and 
to  show  how  this  machinery  works. 

To  maintain  clarity  and  avoid  the  use  of  obfuscating  and  unessential  tech¬ 
nical  details  we  have  restricted  our  reconstruction  results  to  the  classical 
Paley- Wiener  class. 

1.3  Contents  and  notation 

This  paper  is  organized  as  follows: 

Section  2  is  devoted  to  a  brief  summary  of  certain  aspects  of  the  theory 
of  tempered  splines  which  are  germane  to  later  developments  in  this  paper. 
Properties  and  examples  of  Riesz  bases  for  Z,^([— tt,  x])  consisting  of  expo¬ 
nentials  as  well  as  notions  and  results  concerning  mean  periodic 

continuation  are  reviewed  in  Section  3.  The  main  results  together  with  de¬ 
tailed  proofs  are  contained  in  Section  4.  Various  remarks  pertaining  to  the 
material  found  here,  citations,  and  acknowledgements  are  collected  together 
in  Section  -5. 

We  use  standard  mathematical  notation  and  only  alert  the  reader  that 
all  Fourier  transforms  are  interpreted  in  the  distributional  sense  with  the 
normalization  that  when  /  is  an  integrable  function  then  its  Fourier  transform 
/  is  given  by 

m  =  r  /(x)e-"«dx. 

y  — .X> 

As  is  customary,  the  symbols  c  and  C  denote  generic  constants  whose  val¬ 
ues  depend  on  their  occurrence  and  are  usually  independent  of  the  essential 
parameters  in  question. 


2  Tempered  Splines 

2.1  Knot  sequences 

The  type  of  knot  or  sampling  sequences  we  consider  in  this  section  are  se¬ 
quences  of  real  numbers  {rr„}nez  =  <  Xq  <  Xi  <  X2  <  . . .}  which 

satisfy  the  following  properties: 

(i)  A  sequence  {in}  is  sufficiently  dense  if  there  is  a  positive  number  r, 
such  that  intervals  of  length  2r  centered  at  x^.  n  £  Z.  cover  all  of  the 
real  line  IR.  In  other  words 

IR=  U  ^"1  - 

(n)  A  sequence  {in}  is  separated  if  there  i«  a  positive  number  S  such  that 
l^^n  —  >  5  whenever  n  ^  m. 

A  sequence  {in}n€2 's  said  to  satisfy  condition  SDS  whenever  it  satisfies 
both  properties  (i)  and  (ii).  Throughout  this  section  we  will  always  assume 
that  {in}ngz  is  a  sequence  which  satisfies  condition  SDS. 

2.2  Splines  of  even  order 

Splines  were  studied  by  Schoenberg  [23]  in  the  case  of  the  knot  sequence 
{in}  =  Z  and  subsequently  investigated  and  further  developed  by  many 
authors;  see  Section  o  for  more  references.  Here  we  find  it  convenient  to 
restrict  our  attention  to  the  clciss  of  those  splines  which  are  also  tempered 
distributions.  This  is  a  natural  way  of  eliminating  certain  pathological  ex¬ 
amples  and  allows  us  to  use  various  distributional  machinery  including  the 
Fourier  transform. 

Suppose  {in}n6Z  >s  a  sequence  of  real  numbers  which  satisfies  condition 
SDS.  A  tempered  spline  of  order  2k  with  knot  sequence  {i„}ngz  is  a  tempered 
distribution  s  which  satisfies 

(5)  D^*s(i)  =  ^  ar,6(x  -  i„) 

ngz 

d 

where  D  is  the  differentia]  operator  — .  =  DD^~^ .  {a„}„pz  is 

dx  dx^^ 

a  sequence  of  constants,  and  6{x)  is  the  Dirac  measure  at  the  origin.  The 
class  of  all  tempered  splines  of  order  2k  with  knot  sequence  {i„}  is  denoted 
by  SHk{  {.Tn}). 


•  • 


•  •  • 
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The  following  proposition  follows  from  arguments  identical  to  that  used 
to  prove  the  corresponding  fact  in  the  c^ise  {x„}  =  7Z.  Details  may  be  found 
in  Section  2  of  [17], 

Proposition  1  The  following  statements  are  equivalent 

•  s  is  in  SHk{{xn}) 

•  s  is  a  function  which  satisfies  the  following  properties: 

(i)  s  is  in 

(ii)  s  is  of  polynomial  growth.  That  is,  there  are  positive  constants  c 
and  p  such  that 

l-s(x)|  <  c(l  +  |x|)'’. 

{Hi)  On  the  complement  of  {i„),  i?\{x„},  the  function  $  is  infinitely 
differentiable  and  satisfies  D^^s  =  0. 

Thus  S/fjt({x„}„gx)  is  simply  the  familiar  class  of  piecewise  polynomial 
splines  of  order  '2k  which  are  of  polynomial  growth  and  have  the  knot  se¬ 
quence  {Xn}. 

2,3  The  interpolation  problem  and  its  solution 

The  natural  problem  of  interpolation  for  tempered  splines  is  the  following: 

Given  a  sequence  of  real  or  complex  numbers  {yn}n€Z  find  an  element  s 
in  SHk{{xn})  such  that  s{x„)  =  for  all  n  in  Z. 

Since  elements  of  SHk{{xn})  are  of  polynomial  growth  it  is  clear  that 
in  order  for  this  problem  to  have  a  solution  a  necessary  requirement  on  the 
sequence  is  that  it  also  be  of  polynomial  growth.  It  turns  out  that  this 
condition  is  also  sufficient. 

The  following  result  is  a  basic  ingredient  in  the  solution  of  this  problem. 

Theorem  1  (de  Boor  [3])  Suppose  k  is  a  positive  integer  and  {x„}  is  a 
sequence  which  satisfies  condition  SDS-  Then  for  each  integer  m.  m  £  Z. 
there  is  an  element  Xm  in  SHk{{xn})  which  satisfies  the  following  properties: 

1.  For  each  element  in  {i„}  the  function  satisfies 

^rn{Xn}  —  ^m,n 

where  dm.n  h  ronecker  delta. 

•5 
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There  exist  positive  constants  C  and  c  which  depend  only  on  k  and  the 
knot  sequence  such  that 


I 


for  all  m  G  Z. 


Because  of  the  first  property  listed  above  the  functions  are  the  fun¬ 
damental  functions  of  interpolation.  In  analogy  with  classicaJ  polynomial 

1  •  ■  1  1  •  /■  f  1  •  sin7r(x  — m)  ,  ,  , 

interpolation  or  with  the  sinus  cardinalts  basis,  — - ; — .  the  A„  s  are 

ir(i  —  m) 

often  referred  to  as  the  Lagrange  functions  or  the  cardinaJ  functions. 

Consider  any  sequence  {i;„}  which  is  of  polynomial  growth  and  the  cardi¬ 
nal  splines  A^,  m  G  Z,  in  S/fi({x„})  whose  existence  is  guaranteed  by  the 
above  theorem.  The  function 


(6) 


=  Y,  ynA„(x) 


is  well  defined  for  all  real  i  since  the  series  converges  uniformly  on  compact 
subsets  of  the  real  line  IR.  Note  that  s  is  a  member  of  SHk{{xn})  which 
satisfies  s(x„)  =  y„,  for  all  n  in  Z.  Thus  s  is  a  solution  to  the  interpolation 
problem  for  the  data  {pn}-  This  solution  is  unique  since  each  non-zero  spline 
of  order  2k  with  knot  sequence  {x„})  satisfying  s(x„)  =  0  for  all  n  G  ^  must 
have  exponential  growth  in  at  least  one  direction;  see  the  remark  preceding 
Theorem  3  in  [3]. 

We  summarize  these  observations  as  follows; 


Proposition  2  Suppose  k  is  a  positive  integer,  {x„}  is  a  sequence  which 
satisfies  condition  SDS,  and  {(/„}  is  a  sequence  which  for  some  real  number 
p  satisfies 

Pn  =  0(lx„j'’)  OS  n  -»  ±oo. 

Then  there  is  a  unique  element  s  in  SHk{{xn})  which  satisfies 

s{Xn)  =  .Vn 

for  all  n  in  Z.  Furthermore  this  element  s  enjoys  the  representation  (6)  and 
satisfies 

s(x)  =  0(1x1’’)  as  X  -♦  ±00. 


t 
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3  Riesz  bases  and  mean  periodic  continua¬ 
tion 

3.1  Riesz  bases  consisting  of  exponential  functions 

Recall  that  a  basis  of  ^  Hilbert  space  'H  is  called  a  Riesz  basis  if  for 
every  linear  combination  /  =  in  H 

{")  c^|a„|-  <  j|/f  < 

n  n 

where  c  and  C  are  positive  constants  independent  of  /. 

Given  a  sequence  {zn}n€Z  of  complex  numbers  the  collection  of  functions 
defined  by 

(8)  'MO  = 

— cx:  <  (  <  00,  n  £  7Z,  is  denoted  by  Thus 

(9)  ^({-n})  =  {v„}n€Z 

where  the  function  are  defined  by  (8). 

Properties  of  the  collections  5({r„})  ais  bases  where  first  studied  by  Paley 
and  Wiener  [19]  and  subsequently  investigated  by  many  authors;  see  Section 
5  for  more  references  and  further  comments.  Here  our  objective  is  to  provide 
access  to  a  sufficiently  rich  cla,ss  of  examples  of  sequences  {;„}  such  that  the 
corresponding  5({o„})’s  are  Riesz  b£ises  for  tt,  t]).  Toward  this  end  we 
list  two  results:  one  is  formulated  in  terms  of  the  geometric  proximity  of  the 
points  {z„}  to  the  integers  while  the  other  is  formulated  in  terms  of  zeros  } 
of  entire  functions  from  a  special  class  introduced  by  Levin  [15.  16]  which,  in 
analogy  with  the  case  =  'Z.  are  referred  to  as  a’ne  type  functions. 

Theorem  2  (Kadets  [12])  //z„  =  n  +  u;„,  n  £  Z,  where  sup„  |u;„|  <  oc 

and  \Re  u;„|  <  r  <  i  then  is  a  Riesz  basis  for  7r,x]). 

To  state  the  second  result  we  need  to  review  the  following  definition. 

An  entire  function  S{z)  of  exponential  type  is  said  to  be  of  sine  type  if  it 
satisfies  the  following  properties: 

1.  All  its  zeros  are  simple  and  lie  in  a  horizontal  strip.  In  other  words 

{*n}  C  :  \Imz\  <  t/o} 


for  some  positive  number  t/o- 


2.  The  set  of  zeros  {_„}  is  separated,  that  is,  there  is  a  positive  number  S 
such  that  |2m  -  -n|  >  vvhenever  m  ^  n. 

3.  For  some  y  >  jjo  and  positive  constants  m  and  M 

m  <  15(x  +  iy)j  <  M 

for  all  I,  — oc  <  X  <  oc. 

log|5(^y)j 

hmsup  - r-; -  =  TT. 

y—±oa  IJ/I 

Theorem  3  (Levin  [IS],  Golovin  [8])  If  is  the  set  of  zeros  of  a  sine 
type  function  then  the  collection  i’dxn})  is  a  Riesz  basis  for  L^([— jr,  tt]). 

Examples  of  sine  type  functions  may  be  obtained  by  looking  for  repre¬ 
sentations  of  the  form  ^ 

S{z)  =  /  e'^^daU) 

where  a{f)  is  a  function  of  bounded  variation  on  [— z]  having  jumps  at  the 
endpoints  ±z.  If  the  zeros  of  such  a  function  are  separated  then  5(2)  is  a 
sine  type  function.  .A  specific  concrete  example  is  given  by 

S{z)  =  n  sind^-^) 

;=0 

where  the  afs  are  positive  numbers  satisfying  0  <  ao  <  ui  <  .  •  •  <  ay-i  < 

.'V:  its  zero  set  is  L)'’^jQ{N2Z  +  aj)  =  U(..ex{ao  + A'Ar . Q.v-i  choosing 

small  Uj's  shows  that  the  gaps  in  such  a  sequence  can  be  very  large. 

Before  leaving  this  section  it  should  be  mentioned  that  necessary  and  suf¬ 
ficient  conditions  for  the  collection  5({2„})  to  be  a  Riesz  basis  for  t]) 

were  obtained  by  Pavlov  [20j,  see  also  [10].  These  conditions  are  in  terms  of 
zeros  of  certain  entire  functions. 


3.2  Sampling  sequences 

The  kind  of  knot  or  sampling  sequences  we  will  be  interested  in  for  the  rest 
of  this  paper  are  sequences  of  real  numbers  <  xq  <  Ci  < 

X2  <  . . .}  such  that 


(10) 


^({xn}ii6x)  3-  Riesz  basis  for  fd([-T,T]). 


~ry 
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VV'e  call  such  sequences  RRB  sequences  or  say  that  it  satisfies  condition  RRB. 
Thus.  {  jT  is  an  RRB  sequence  if  it  is  a  sequence  of  real  numbers  and  the 
corresponding  sequence  of  exponentials  enjoys  property  (10). 

It  is  clear  that  every  RRB  sequences  must  be  separated,  that  is.  it  must 
have  property  (ii)  in  Section  2.1.  Otherwise  it  is  easy  to  check  that  the  left 
hand  side  inequality  of  (7)  fails. 

Furthermore  it  is  a  direct  consequence  one  of  Beurling’s  celebrated  results, 
see  [2.  14],  that  every  RRB  sequence  {ar„}„gz  must  satisfy 

(11)  limintf  int  , 

r— oo  — •x><a<3o  p  j 

where  (a.  a  +  r)  denotes  the  interval  a  <  i  <  a  +  r  and  denotes  the 
number  ol  elements  in  the  set  Q.  It  is  clear  that  any  sequence  which  satisfies 
(11)  must  be  sufficiently  dense,  that  is,  it  must  satisfy  property  (f)  in  Section 
2.1.  Otherwise  it  is  easy  to  check  that  th'*  right  hand  side  of  inequality  (11) 
is  zero. 

These  observations  are  summarized  as  follows:  Every  RRB  sequence 
satisfies  condition  SDS. 


A 


3.3  The  operator  of  mean  periodic  continuation 

Suppose  that  {jr„}  is  an  RRB  sequence.  Then  each  function  h  in  Z,^(i-~.  :r]) 
enjoys  the  representation 

(12)  /*(0  =  E  «n  -IT  <  ^  <  TT, 

where  the  sequence  of  coefficients  {a„}  satisfies 

(13)  c||{a„}||<qx)  <  ||/8||iJ((_,.;r|)  <  C||{an}||rqz) 

where  c  and  C  are  positive  constants  independent  of  h.  Note  that  the  series 
on  the  right  hand  side  of  (12)  defines  a  prolongation  of  the  function  k  to  the 
whole  real  line.  We  call  this  prolongation  //,  Thus 

(M)  //(O  =  E  '*’>  <  (f  <  oc, 

where  the  sequence  of  coefficients  {on}  is  the  same  as  that  in  (12).  By  virtue 
of  f  1.'))  the  series  f  1  1)  is  locally  convergent  and  the  function  is  locally 
in  . 


q 


•  ••••••• 
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Consider  the  linear  operator  A  :  k  Ah  which  maps  I"([— tt,  t])  into 
itself  and  is  defined  via 

(15)  Ah{0  =  E  ««  -Jr  <  ^  <  x, 

where  the  coefficient  sequence  {a„}  is  uniquely  determined  by  h  via  (12). 
Because  the  functions 

(16)  <p„(0  =  n  G  S 

are  a  Riesz  basis  for  L^([-7r,  tt])  the  operator  A  is  well  defined  and  invertible 
on  L^([— TT,  tt]).  Similarly  the  integers  powers  of  .4,  A-'  =  A.A-’"\  j  E  Z. 
may  be  defined  via 

(17)  A^h(0  =  E  -t  <  ^  < 

ne2 

Again,  because  the  collection  of  functions  £’({i„})  =  defined  by  (16)  is 
a  Riesz  basis  the  operator  norms  of  A-'  are  uniformly  bounded,  namely, 

(18)  ||A-'/i||£2((_,r,T|)  <  C'Plk2((_;r,r))-  j  E  Z 

where  C  is  a  constant  independent  of  j  and  h. 

Note  that  A-'/i  can  be  e.xpressed  in  terms  of  the  prolongation  H  as  follows: 

A^h{f,)  =  H((,  +  'lltj),  -TT  <^  <  TT. 

For  this  reason  A  is  referred  to  cis  the  prolongation  operator  ov  the  operator 
of  mean  periodic  continuation  with  respect  to  the  collection  5({xn}). 

We  will  also  need  a  representation  of  the  adjoint  .4“  of  .4.  To  this  end 
consider  the  Riesz  basis  which  is  dual  to  in  /.’([-t.t]).  Then  if 
h  in  Z;^([-7r,  tt])  enjoys  '.he  representation 

(19)  h{a=Z<^n(0- 

where  the  coefficient  sequence  {a'}  is  uniquely  determined  by  h,  it  is  clear 
that 

(20)  A-hio  =  E  <  ^ 

and.  more  generally,  for  all  j  G  Z.  (A'Y  =  (.A-')"  =  .-A'-'  has  the  representa¬ 
tion 

(21)  A-^h{f)  =  E  <  <  ^  <  ^- 

Since  it  will  be  u.sed  often  in  what  follows,  we  remind  the  reader  that  the 
operator  norms  of  {.4"^}  are  also  uniformly  bounded. 


•  •  • 
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4  The  recovery  of  band-limited  functions 

4.1  The  Paley- Wiener  class  and  spline  interpolation 

If  0  <  <  oc  the  Paley-Wiener  class  PWp  consists  of  those  functions  in 

i^(-oo,oo)  whose  Fourier  transforms  have  support  in  the  interval 
In  other  words 

FIVj  =  |/ e  F^(-oc,  oc)  :  supp /  C  [-;3,  aJ]|. 

We  remind  the  reader  that  such  functions  /  are  restrictions  to  the  real 
line  of  entire  functions  F{z)  of  exponential  type  /3.  In  particular  such  f's  are 
continuous  and  well  defined  point  wise. 

Suppose  {x„}  is  a  sequence  which  satisfies  condition  RRB.  If  /  in  PW,^ 
then  the  sequence  of  values  is  in  Define  5;-/  to  be  the  spline 

in  5^fjt({a;„})  which  interpolates  /  on  {in},  namely, 

Skfi^n)  =  fi^n)  for  all  n  e  Z. 

In  view  of  Proposition  2  the  function  5{./  is  well  defined;  that  is.  Skf  exists 
and  is  unique. 

The  remainder  of  this  paper  is  devoted  to  showing  that 

lim  Skf(i)  =  fix) 

uniformly  and  in  T^(-oo,oc).  Thus  in  what  follows  we  always  assume  that 
{in}  is  an  RRB  sequence,  /  is  in  PW^,  and  S^f  is  the  element  in  SHki{xn}) 
which  interpolates  /  on  {in}- 

4.2  Elementary  properties  of  Skf 

In  view  of  the  e.xponential  decay  of  the  An’s  and  the  fact  that  5({i„})  is  a 
Riesz  basis  for  Z/^([— i,  i])  it  should  be  clear  that  the  mapping  /  — ,  5^/  is  a 
linear  transformation  from  PW^  to  Z,^(-oc.ool  which  satisfies 

(22)  <  c||{/(i„)}l|,.(^,  <  C\\fU^^.^,^^ 

where  c  and  C  are  positive  constants  independent  of  /. 

Since  D^'‘Skf(x)  =  a^Six  —  i„)  we  have 

e'^sTfif)  = 


(23) 


where 

(24) 


nex 

SkfiO  is  in  Z,^(-oc,oo)  implies  that  is  locally  in  L'  and,  in  par¬ 

ticular,  is  in  L^([— TT,  tt]).  Because  of  (24)  and  the  fact  that  {a:„}  is  an  RRB 
sequence  we  may  conclude  that 

where  C  is  a  constant  independent  of  j  and  that  <  oo-  For  later 

reference  we  denote  the  restriction  of  to  the  interval  [— x,  x]  by  in 
other  words  ijjk  =  ^fc|[_,r,)r]  and  if  {2j  —  l)x  <  ^  <  {2j  +  l)x  then 

(25)  -  2tj) 


for  all  j  in  Z  where  is  the  prolongation  operator  defined  in  the  previous 
section. 

In  view  of  (23)  we  may  write 

(26)  ^ 


II  W(Olli.,-,o.oc)  =  +  E  /■' 

which  by  virtue  of  (25)  may  be  re-expressed  as 


(2;+i)^|«l;^(^)|2 


(2;-l)fr  I 


(27)  ||•5'^/(^)|li2,_.„,^)  =  l|5')t/(Ollf,2([_,,,j)  +  E  /  |/f  , 

;e*\{0}'^-’'  'I** 

Proposition  3  The  function  xpk  satisfies  the  relation 

for  —IT  <  ^  <  X  where  ,4  is  the  prolongation  operator  and  .4*  is  its  adjoint. 


Proof:  Write 

r  m)e‘^''^d^  =  2x/(xJ  =  2x5V(x„) 

J-^K 


and  hence 


=  /“ 


'  /:{s-((?^)<«}'“^ 


(29)  £  /(«e-<i«  .  £  {  E  ((JTW(''‘'«’‘)('»)K)}'-*‘'f' 
Since  (29)  is  true  for  aJl  n  in  ^  and  {-Cn}  is  an  RRB  sequence  it  follows  that 
/(fl  =  E/-'(7rTi-p-«*(,))(£) 


which  is  the  desired  result. 


For  later  reference  we  re-express  (28)  as 

(30)  iV(0+  E  +  = 

;6Z\{0}  ^  ^  ' 

for  -ff  If  we  use  the  notation  s*  to  denote  the  restriction  of  Skf 

to  the  interval  [— t.  jt],  in  other  words  =  5t/|[_r,;r],  then  (30)  may  be 
expressed  as 

(31)  5i(0+  E  '4'^((^-2;y“(A4(-)''s,(c.-)l)(r,))(0  =  /(a 

jez\{o}  ^  ^ 

for  —  ;r  <  <  TT.  This  expression  may  be  abbreviated  to 

(32)  {I  +  Bk:\h)h  =  f 


where  ih  and  /  are  elonients  of  //■*([— r,  jt])  and  /,  Bh,  Mk  are  linear  operators 
on  Z,*([-7r,  ir])  defined  by 

BkgiO  =  Z  + 

jex\{o)  ^  ^  ' 

A/t5(o  =  {a^?‘‘9{a 

and  I  is  the  usual  identity  operator  lg{0  =  9i0-  Note  that 

l|5fc5lU>((-^.^l)  -  {  5]  (2|i|  - 


'  i€  A{o} 


or,  more  simply, 

(33)  ||Bt5||lJ((_,,,))  <  Cll5||L»([-»,rl) 

where  C  is  a  constant  independent  of  k  and  g. 

Lemma  1  If  f  is  in  PW^  and  0*  is  related  to  f  via  (2S)  then 

(S4)  ll’/’fclU^d-T.irl)  <  ’r*‘ll/lk’(l-ir.)r]) 

Proof:  Take  the  Z/^([— t.t])  scalar  product  of  the  function  0;.  with  both 

sides  of  (28).  This  results  in 

Since  the  summands  on  the  left  hand  side  of  (35)  are  non-negative  we  have 
(36)  (^,  00  <  |(/,  0O|. 

Inequality  (36)  together  with 


|(/,  001  ^  ll/lU>(I-)r,rl)||0fc||/:»((-»r.itI) 


imply  that 


^  ll/i!/i>((-<r,)rl)110fclU>((-ir./r]) 
which,  upon  simplification,  is  the  desired  result. 


Lemma  2  If  f  is  in  PW^  then 

ll5i/llLJ([-ff,r))  <  C'll/l|L2([-)r,T]) 
where  C  is  a  positive  constant  independent  of  f  and  k. 
Proof:  By  virtue  of  (30)  we  may  write 

l|-5fc/lli3([_;r,T])  <  ll/lU^d-T.,]) 

(38)  +  E 

jez\{o}  ^  ^  ^  ^ 


("2|j|  -  l)2i  ■  ^ 


where  C  is  a  constant  independent  of  j,  fc,  and  /.  In  view  of  (34)  the  last 
inequality  together  with  (38)  imply  the  desired  result.  ■ 


Proposition  4  The  mapping  f  Skf  from  PW.„  to  L^([— oc.oo])  is  bounded 
uniformly  with  respect  to  k.  In  other  words,  if  f  is  in  PW,.  then 

(39)  l|•!>t/||z,2(-0O,0o)  <  C'II/IU2(-«),00) 
where  C  is  a  positive  constant  independent  of  k  and  f. 

Proof:  In  view  of  Plancherel’s  formula  (39)  is  equivalent  to 

(40)  \\^f\\m-oo.oo)  < 

To  see  (40)  use  (27),  the  identity 


and  the  fact  that 


to  write 


(^1)  +  (2|;i  _ 


J€Z\{0) 


I 


Sinc'3  the  constant  C  in  (41)  is  independent  of  /  and  k,  inequalities  (34)  and 
(37)  together  with  (41)  imply  the  desired  result.  ■ 


4.3  Behavior  of  5t/  as  k  —*  oc.  ' 

Theorem  4  If  f  is  PW^  then 

(4^)  lim  11/  -  =  0 

«— *>00 

Proof;  By  virtue  of  Plancherel’s  formula  (42)  is  equivalent  to 
143)  ^lim  11/  -  5i/||i,2(_^,«>)  =  0. 

To  see  (43)  use  the  fact  that  /  vanishes  outside  the  interval  [— tt,  t]  and  write  l 

11/ ~ -^fc/ Hi’ (-»,»)  =  \\f  ~  ^kf\\h{i-^,K]) 

(44) 

+  H  l|•^»:/lli’(((2}-l)^^,  (22  +  1)t))- 

jez\{0} 

Estimate  the  size  of  the  first  term  on  the  right  hand  side  of  (44)  as  follows:  * 

Let  Sit  denote  the  restriction  of  St/  to  the  interval  [— tt,  tt]  and  recall  relation 
(32),  which  is 

(/  +  BkMk)sk  =  /. 

In  view  of  Proposition  4  it  follows  that  the  operators  /  +  BkMk  are  invert-  t 

ible  as  mappings  from  L^([— n-,  rj)  to  itself  and  the  inverses  are  uniformly 
bounded,  namely 

(45)  ||(/  4-  5ltMt)"'|(t2(j_„,,r))-£,’([-,r,^j)  <  C 

where  C  is  a  constant  independent  of  k.  Thus  we  may  write  * 

/  -  it  =  /  -  (  /  +  BkMkr'f  =  (  /  +  BkMk)-^BkMkf 

or  more  directly 

(46)  f-sk  =  {l  +  BkMk)-^BkMkf.  ► 

Recall  that  the  operator  norm  of  Bk  is  bounded  independent  of  k,  see  (33). 

This  together  with  (45)  and  (46)  imply  that 

(47)  11/  -  i’t/llt’d-T.r))  <  C'l|-V|t/||i,2([_c.,T)) 


t 


•  •  •  • 
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which  is  the  desired  estimate  for  the  first  term. 

Estimate  the  size  of  the  remaining  terms  on  the  right  hand  side  of  (44) 
a.s  follows:  Use  reasoning  similar  to  that  used  to  obtain  (27)  to  write 


(48) 


-l)>r,(2j+l)ir]) 

J€Z\{0} 


S  f 


je^fo} 


U  + 


d^. 


Note  that 
(49) 

where,  as  above,  is  the  restriction  of  5*/  to  the  interval  [— :r,  x]  and 
MkSkiO  =  (^/’r)^*4(0-  Thus 


(50) 


U  +  2^jP‘‘ 


A^MkhU) 


and  the  term  on  the  right  hand  side  of  (48)  is  equal  to 


which  is  dominated  by 


{ (21;"-  • 


Hence  we  may  write 

(51)  l|•5£/||t^(((2,-lK,  (2j+l)ir))  ^ 

2€Z\{0} 

where  C  is  independent  of  /  and  k.  Now 

(52)  l|445tllL2((-T,T))  <  ~  /)|U2([-r,!r)) 


and  clearly 

||T/fc(^*:  -  /)ltt2((-r,rl)  <  ll^fc  -  f]\m[-n.r\)- 
In  view  of  (47)  the  last  inequality  may  be  replaced  by 

(53)  ||Mt(.9r.  -  /)||z.2((-r.r])  <  (^\\Afkf\\m{-n.K]) 

where  C  is  a  constant  independent  of  /  and  k.  Combining  (51).  (52),  and 
(53)  results  in  the  estimates 

(•54)  \\-^h^k\\L'H{-,r.K\)  <  ^ll•''^4/||z,*((-l^.^]) 


> 


I 


» 


) 


I 


•  •  •  • 
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and 

(■5^)  ll•^fc/llL2([(2J+l)lr,  (2>+l))r])  - 

J6^{0} 

where  C  is  a  constcint  independent  of  /  and  k. 

Estimates  (47)  and  (55)  when  substituted  into  (44)  lead  to 

(56)  11/  -  SkfWii^-oo.oo)  < 

where  C  is  a  constant  independent  of  /  and  k.  Inequality  (56)  implies  the 
desired  result  since  Mkf{^)  =  (//7r)^*/(^)  and  the  Lebesgue  dominated  con¬ 
vergence  theorem  imply  that 

lim  =  0. 

fc— OO 


Inequality  (56)  is  interesting  and  worth  re-statement  as  a  corollary  of  the 
above  argument. 

Lemma  3  If  f  is  in  PW^  then 

(57)  11/  -  5,/lli.(_^.«)  <  C7  £  (I  j  ‘  m  di 

where  C  is  a  constant  independent  of  k  and  f. 

The  technique  used  in  the  proof  of  the  theorem  leads  to  other  convergence 
results.  For  e.xample.  to  estimate  the  pointwise  difference  between  /  and  Skf 
we  may  write 

Skfix)  -  fix) 

(58)  =  ±1  /  ^m)e‘^^dA 

ItC  L  J-ft  J 

The  first  integral  on  the  right  hand  side  of  (58)  can  be  estimated  via  the 
Schwartz  inequality  and  (47)  to  get 

(59)  I  r  (STm)  -  /(/))e-"«d^  < 

I  7  — TT 

where  C  is  a  constant  independent  of  /  and  k.  To  estimate  the  second 
integral  on  the  right  hand  side  of  (58)  use  (26)  and  rea.soning  similar  to  that 
used  to  get  (27)  and  write 


/.  SJ{ae-'^d^=  ^  T- 


;€Z\{o}  ■ 


A^U'kif)  ix{(+2r:)ie 
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Now  use  (49)  and  (50)  to  see  that  if  /,  is  the  j-th  term  on  the  right  hand 
side  of  (60)  it  can  be  expressed  £is 

fj  =  /_"  -  2j) 

Using  the  Schwartz  inequality  results  in 

(61)  1/,|  <  l|AU¥fcS,(0||L^((-.,M)  •  {  /_'  If  '  2j 


Since  the  first  and  second  terms  in  the  product  on  the  right  hand  side  of  (61) 
are  dominated  by 

\/2w 

C\\Mkf\\m[-n,r])  and  _Y)-2k 

respectively  (the  first  follows  from  (54)  and  the  second  by  direct  estimation) 
we  have 

(®2)  |/ji  <  (2|j|  - 

where  C  is  independent  of  j.  k.  and  /.  Thus 


j€X\{0> 


r.^1) 


■  jez\{o} '-IjI 

which  may  be  simplified  to 

(63)  I  f  <  C||.V/J||tq[_,.,), 

where  C  is  a  constant  independent  of  k  and  /. 

Combining  (59)  and  (63)  to  estimate  the  integrals  on  the  right  hand  side 
of  (58)  leads  to  a  bound  on  |5;t/(z)  —  /(z)|.  We  summarize  this  as  follows. 

Theorem  5  //  /  is  in  PW^  then  for  all  .r  in  IR 


(64) 


\S,f(.r)-  f{x)\<c\^£^  no 


whtre  C  is  a  constant  independent  of  k  and  f .  It  follows  that 
(65) 

uniformly  on  [R. 


lim  Skfix)  =  ffr) 

fe— OO 
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Estimates  (57)  and  (64)  lead  to  even  better  convergence  results  if  the 
support  of  /  is  properly  contained  in  the  interval  (  — tt,  tt).  The  rate  of  con¬ 
vergence  is  geometric  and  depends  on  the  size  of  the  support.  In  particular 
if  the  support  of  /  is  contained  in  the  interval  [— /?,  0]  where  0  <  <  tt  then 

the  integral  on  the  right  hand  side  of  (57)  and  (64)  is  dominated  by 

(f)'  11/11.^,-00,00)  =  ll/IU^(-.oo.oo). 

We  summarize  this  remark  cis  follows: 

Corollary  1  If  f  is  in  PWj}  for  some  0  which  satisfies  0  <  0  <  tt  then 

(6^)  11/  -  •^i/llz;^(-00,00)  <  C"  ||/||.2(_oo,2o) 

for  all  p  which  satisfy  2  <  p  <  oc  where  C  is  a  constant  independent  of  k 
and  f. 


5  Remarks  and  acknowledgements 

In  the  next  three  paragraphs  we  collect  several  remarks  and  indicate  further 
references  which  may  be  helpful.  No  attempt  is  made  to  be  exhaustive. 

The  theory  of  splines  and  their  applications  has  undergone  extensive  de¬ 
velopment  in  the  past  thirty  or  so  years.  In  addition  to  the  references  cited 
above  we  mention  that  [25]  is  a  relatively  recent  treatise  on  the  subject  which 
contains  an  extensive  list  of  references.  We  note  that  our  version  of  Propo¬ 
sition  2  is  an  ecisy  consequence  of  the  material  in  [3]  and  bring  attention  to 
the  fact  that  the  statement  of  a  version  with  a  less  restrictive  condition  on 
the  sequence  of  knots  {i„}  can  be  found  in  [11]. 

A  history  and  contemporary  exposition  of  the  theory  of  Riesz  bases  con¬ 
sisting  of  exponentials,  including  the  results  quoted  in  Section  3.1,  is  given  in 
the  survey  [lOj.  See  also  the  textbook  [27],  More  recent  surveys  and  lists  of 
references  may  be  found  in  [I.  6].  Material  on  mean  periodic  functions  may¬ 
be  found  in  [13,  26]. 

The  fact  that  various  classes  of  band-limited  functions,  or  entire  functions 
of  exponential  type,  can  be  recovered  from  their  samples  on  certain  discrete 
sets  is  well  known.  The  general  subject  area  is  often  referred  to  as  sampling 


I 


•  •  • 


•  •  •  • 
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theory.  For  surveys  which  include  the  result  involving  condition  (11)  men¬ 
tioned  in  Subsection  3.2  and  a  list  of  further  references  see  [1.  6];  indeed  most  • 

of  the  work  cited  in  the  previous  paragraph  is  directly  related  to  this  sub¬ 
ject.  Interesting  surveys  which  include  various  e.xtensions.  generalizations, 
error  estimates,  numerical  methods,  and  further  references  can  be  found  in 
[5,  7,  9]. 

This  work  was  completed  during  the  spring  semester  of  1993  while  * 

Lyubarskii  enjoyed  a  visiting  position  with  the  Department  of  Mathematics 
at  the  University  of  Connecticut;  he  is  grateful  to  the  Department  for  this 
opportunity.  Madych  was  partially  supported  by  D.4RP.A  Grant  AFOSR-90- 
0311. 

t 
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Orthogonality  Criteria  for  Compactly 
Supported  Scaling  Functions 

Karlheinz  Grochenig  * 


Abstract 

V\e  give  a  new  simple  proof  of  several  criteria  for  the  orthogonality 
of  the  integer  translates  of  a  scaling  function. 

-A.  scaling  function  is  a  function  €  V^iJR)  that  satisfies 

(1)  c*  g  -  i) 

kez 

for  some  integer  q,  |?1  >  1,  and  coefficients  c*  G  €,  and  such  that  the  integer 
translates  —  k),k  6  Z,  are  an  orthonormal  system.  In  order  to  construct 
a  compaclty  supported  scaling  function  >p  [5],  one  chooses  a  finite  sequence  of 
non-zero  coefficients  c*  €  (T,  c*  =  0  for  |il'|  >  N,  such  that 

(-)  ?  =  6t  for  all  ££  Z 

kez 

and 

(3)  Z  =  1 

kez 

If  m(^)  denotes  the  corresponding  trigonometric  polynomial 

(■1)  m{0  =  Z  > 

k=-N 

properties  (2)  and  (3)  are  equivalent  to 


(5) 


for  all  ^  ^  IR 


and 

(6)  ^”(0)  =  1  . 

*The  author  acknowledges  the  partial  support  by  a  DARPA  Grant  AFOSR-90-0.311. 
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The  scaling  function  ^  is  then  obtained  from  the  infinite  product 

(7)  ‘PiO  =  n 

jr=l 


and  Fourier  transform  <p(x)  =  //j  The  infinite  product  converges 

uniformly  on  compact  sets  to  (p  €  has  compact  support,  and 

satisfies  the  scaling  relation  (1),  see  for  instance  [5],  Ch.  6. 

The  question  whether  the  integer  translates  of  such  a  ip  are  mutually 
orthogonal  is  more  subtle.  Conditions  (2)  and  (3)  are  necessary,  but  not 
sufficient  to  imply  orthogonality. 

In  this  article  we  give  a  new  and  simple  proof  of  several  orthogonality 
criteria  obtained  by  A.  Cohen,  W.  Lawton,  J.  P.  Conze  and  A.  Raugi,  and 
Q.  Sun. 

Since  “almost  all”  sequences  satisfying  (2)  and  (3)  lead  to  orthogonal 
translates,  we  characterize  the  case  of  non-orihogonal  translates. 

The  non-tri  vial  part  of  the  proof  will  consist  in  the  analysis  of  the  operator 


By  (5)  the  constant  1  is  always  an  eigenfunction  for  the  eigenvalue  1. 


Theorem  1  Let  m{^)  be  a  trigonometric  polynomial  satisfying  (5)  and  (6) 
and  p(^)  =  njii  Then  the  following  statements  are  equiv¬ 

alent: 

(A )  The  translates  p(z  —  k),  k  E  Z,  are  not  orthogonal. 

(B)  There  exists  a  nonconstanl  nonnegative  trigonometric  polynomial  f 
so  that  Tf  =  f. 

(C)  There  exists  a  ^  E  (0,  J)  so  that  =  ^[mod  1)  for  some  N  >  1  and 

=  1  /o'"  o/f  j  >  0. 


(D)  There  exist  integers  N,s  >  0,  I  <  s  <  q'"^  —  i  so  that 

(  q^s 

- T  ^  q)  ~  ^  ^  0  and  £  =  1,2, . . .  ,q  —  I 

(E)  There  exists  q  E  (1,2)  with  the  following  property:  For  all  k  E  Z 
there  is  a  j{k)  >  1  such  that 


I 


I 


I 


» 


» 


I 


» 


I 


> 
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Remarks  1.  The  equivalerce  (A)  (E),  due  to  A.  Cohen  [l],  was 

the  first  known  characterization  of  orthogonality.  It  is  usually  formulated 
positively  in  the  following  way;  The  integer  translates  of  -p  are  orthonormal  if 
and  only  if  there  exists  a  compact  set  K  C  IR  which  contains  a  neighborhood 
of  0  and  satisfies  Utez(^  +  K)  =  IR,  (A:  +  if)  fl  =  0  for  A;  /  0,  such  that 
0  for  all  j  >  1,  and  ^  €  if  (“m  satisfies  Cohen’s  condition”). 

The  characterization  (.\)  (B)  was  found  independently  by  W.  Lawton 

[8]  and  J.  P.  Conze  and  Raugi  [4].  The  arithmetic  characterizations  (C) 
and  (D)  appear  in  A.  Cohen  and  Q.  Sun  [2],  and  (.A)  (C)  can  also  be 

found  in  [4].  The  analysis  of  eigenfunctions  of  T  also  figures  prominently  in 
[3]. 

2.  The  equivalence  (,A)  <=>•  (B)  <=>  (E)  e.xtends  to  more  general  filters  m 
and  also  holds  in  higher  dimensions  [i,  8,  9]. 

3.  Our  contribution  is  a  new  organization  of  the  proof  which  makes  all 

but  one  implication  trivial  and  a  new  argument  for  the  crucial  step  from  (B) 
to  (C).  Compared  with  the  proofs  in  the  references  above  or  in  [5]  our  proof 
is  simpler  and  much  shorter.  ■ 


Proof:  We  first  remark  thai  by  Poisson’s  summation  formula 

(8)  =  Z  (  L  . 

k€Z  k^Z  ' 

If  p  has  compact  support,  then  /(^)  Z*ez  \‘p{^  ~  's  a  trigonometric 
polynomial,  and  the  translates  p[x  —  k)  are  orthonormal  if  and  only  if  /(^)  = 
1  for  all 

(.A)  (B).  If  {p{x  ~  k),  k  G  2Z}  is  not  orthonormal,  then  /(^)  is  a 

nonconstant  nonnegative  trigonometric  polynomial.  Using '^(^)  =  m{^)p{^), 
we  obtain 


(9)  /(O 


=  =  Z  Z 

k  i=0  t^Z  ^ 


ZJ.'S 

) 

(B)  (C).  Set  m„(f)  =  n"Co  Then 

00. 


*=0 


as  is  easily  verified  by  induction.  Since  Tl  =  1,  we  also  obtain 


(11) 


=  1  for  all 


3 


I 


Now  assume  that  T f  =  f  for  some  nonconstant  nonnegative  trigonomet¬ 
ric  polynomial.  Then  one  of  the  extrema  is  assumed  at  an  interior  point 
foiO  <  Co  <  1-  Without  loss  of  generality  we  may  assume  /(^)  <  /(Co)  for 
all  f  €  fR.  The  hypothesis  T"/  =  /  for  all  n  and  (11)  yield 

(.2)  0  =  /  -  r-/  =  |m(^)|’(/K„)  -  /(^)) 

As  all  terms  in  (12)  are  nonnegative,  either  /(Co)  =  /(^^)  or  rn„(^'^)  =  0 
holds  for  each  k. 

Let  M  be  the  finite  set  of  maxima  of  /  and  C  =  card  M.  Then 


(13)  {ry  :  7?  =  and  /  0}  C  tW 

Consequently  for  each  n  the  sum  in  (11)  contains  at  most  C  non-zero 
terms  and  there  is  a  0  <  k„  <  q",  such  that 


(14) 


m„ 


( 


Co  +  \  ^ 

q"  J  ~ 


—  t 


We  can  now  choose  an  infinite  subsequence  n(j)  so  that 


(15)  -^;^  =  ^€.W,C^0,1  for  all  t. 

Since  ^'’''^C  =  Co  +  ^n{.)  =  Co(  mod  1),  there  exists  an  N  >  1  such  that 
q^^  =  mod  1),  e.g.  j\'  =  n{2)  —  n(l). 

Finally,  since  =  q-’^{  mod  l)  for  all  l,j  >0  and  since  the  sequence 

|m„(C)|,n  =  1,2,...,  is  decreasing,  we  obtain  for  L  >  1 


(16) 


N-l  .iA'-l  ^i/L 

n  lm(?-'OI  =  (  n  |rn(?^C)l)  > 

j=0  j=0  ' 


>  lim-o 


n(i)-l 


l/L 


n  |m(?-'C)  =  liimv-oc  m„(.)(0 

j=o  '  ' 


l/L 


>  C”^ 


.As  L  was  arbitrary,  we  obtain  the  desired  conclusion  |m(g-’^)|  =  1  for  all 
/  >  0. 

(C)  =>  (D).  .As  q^^  =  ^  (mod  1),  we  have  q"^  =  ^  +  s  hr  some  s,  1  < 

s  <  q^  —  1  and  thus  ^  By  (5)  and  by  the  assumption  |m(?-'C)|  =  1 

we  obtain  -f.  ^)  =  0  for  ^  =  1, ...  9  —  1. 

(D)  =»  (E).  Set  jy  =  iphn  +  verify  that  for  each  k  ^  Z  there 

is  a  J  £  Z.  J  >  0,  such  that 


(17) 


'7  +  ^  ^  , 

qJ+i  ,.v  _  1 
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for  some  r€  1},  /  €  —  1}  and  M  £  Z.  If  (17)  holds, 

then  by  hypothesis  =  0. 

Given  k.  €  Z ^  we  consider  the  sequence 


1  ■  1 

aj  =  q-^(-s'£.q'‘'^  +  k+\ 


t  +  i 

^  UN  1  \  j  '  7  1 )  -  *  * 

q^\q  -  1)  9-' 


where  [x]  denotes  the  smallest  integer  <  x.  The  accumulation  points  of  {uj} 
are  the  numbers  —  -% ^  Z,  r  —  0. 1, . . . ,  jV  —  1.  Therefore  there  is  a  largest 
integer  J  >  0,  such  that  aj  €  Z,  but  uj+i  ^  Z,  i.e.,  aj+i  =  '-  +  M  for  some 
/€{l,2,  —  1}  and  M  €  Z.  Rewriting  oj+i  as 


q->*^  V 


^[±]+N-J. 
q^  -  1 


yields  (17). 

(E)  =>  (A).  If  for  all  i  +  ^))  =  0  for  some  j  >  1,  then 

<f{q  +  k)  =  n^i  'n{q~^{q  +  i))  =  0  for  all  k.  Thus 


Y.  =  0 

kex 


and  the  translates  .p{x  —  n)  cannot  be  orthonormal  by  the  remark  at  the 
beginning  of  the  proof.  ■ 


The  theorem  furnishes  a  quick  proof  for  the  classification  of  self-similar 
tiles  and  Haar  bases  of  L^{,IR)  which  was  obtained  in  [6).  Given  q  €  Z,  \q\  > 
1,  let  V  =  (Aiq,  . . . ,  be  a  complete  residue  system  modulo  q,  i.e. 

ki  =  I  (mod  q),  and  set  c*  =  1/g  for  fc  €  X>  and  C),  =  0  otherwise.  Set 

(18)  Q  =  Q{V)=  {xeR:x  =  Yq-’e„i,&  D}  , 

j=i 

then  Q  is  a  compact  and  self-similar  set  satisfying 

qQ  =  uV<  +  Q)  ■ 

1=0 

Therefore  its  characteristic  function  <p  =  XQ  '•he  (unique)  solution  of  e>(i)  = 
IIITg  <i>{qx  —  ki).  {(i)(x  —  k),  k  £  Z]  is  orthonormal,  if  and  only  if  Q  n(i  Q) 
has  measure  0  for  all  fc  ^  0,  in  other  words,  if  and  only  if  Q  tiles  R. 

In  the  following  gcd  denotes  the  greatest  common  divisor  of  a  given  set 
of  integers. 


Theorem  2  ([6]  )  The  integer  translates  k  +  Q,k  €  Z,  are  mutually  dis¬ 
joint  (up  to  sets  of  measure  zero)  if  and  only  if  gcd  i,j{ki  —  ij)  =  1. 

Proof:  By  the  orthogonality  criterion  /  Xe(^)xo(^  —  k)dx  ^  0  for  some 

i  ^  0,  if  and  only  if  there  exists  ^  €  (0, 1),  such  that  =  ]  for  all 

i  >  0.  In  this  case 

1^(01 = I- r  = 1 

?  i=0 

implies  that  =  e^"'“  for  some  a  €  [0,1]  and  /  =  0, 1, —  1 .  Since 

kt^  =  a  (mod  1),  it  follows  that  (i,-  —  kj)^  €  Z  for  all  i,j  =  0,\, . . .  ,q  —  1 
and  by  taking  appropriate  linear  combinations  also  gcd  —  kj)^  €  Z. 
Since  0  <  ^  <  1,  we  obtain  gcd  ,j{ki  —  kj)  >  1. 

Conversely,  if  d  =  gcd  —kf)  >  1,  then  lrn(9-’^)P  =  1  for  all  j  >  0.  ■ 

It  is  now  easy  to  obtain  explicit  compactly  supported  wavelet  bases  that 
are  analogous  to  the  ordinary  Haar  bases.  See  [7]  for  more  details. 

CoroUary  1  ([7])  Let  q,ki,Q  be  as  above  and  assume  that  gcd  —  kj)  = 

1.  Let  U  =  . 9-i  be  a  unitary  q  x  q-matrix,  such  that  Uoj  = 

;  =  0, . . . ,  ^  —  1.  Define 

9-1 

(19)  rvfix)  =  Y.  -  ^i)  fori=  \,...,q-\  . 

j=o 

Then  the  collection  of  functions 

q-’l^iffiq^x  —  t)  j,  k  E  Z,  i  =  \, - 9  —  i 

is  a  complete  orthonormal  basis  for  L^{BT). 
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